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ABSTRACT 

Nonlinear finite element analysis was applied to various reinforced concrete beams using a set of constitutive models 
established in the modified fixed-angle softened-truss model (MFA-STM). The model was implemented by modifying the 
general-purpose program FEAPpv. The model can take account of the six important characteristics of cracked reinforced 
concrete: (1) the softening effect of concrete in tension-compression; (2) the tension-stiffening effect of concrete in tension; 
(3) the average stress-strain curve of steel bars embedded in concrete; (4) the shear modulus of concrete; (5) the aggregate 
interlock; and (6) dowel action. The comparison shows the aggregate interlock and dowel action can reduce the 
overestimation of the shear capacity of high strength reinforced beam, especially the high strength reinforced deep beam 
without web reinforcement. Moreover, the model is suitable for being implemented numerical procedures due its simplicity. 

KEYWORDS: high-strength concrete; shear; aggregate interlock; dowel action; finite element. 



INTRODUCTION 



DESCRIPTION OF MODEL 



The nonlinear finite element method has developed 
into an important tool for the analysis of the complex 
concrete structures. This technique is very helpful to 
understand the formation and propagation of cracks and 
the mechanism and process of failure. Future development 
of the nonlinear finite element method lies primarily in 
the improvements of the constitutive models of materials. 
Two behavioral models were developed for the analysis 
of concrete structures subjected to shear: the rotating- 
angle softened truss model (RA-STM) 13 and the fixed- 
angle softened truss model (FA-STM). 3-6 The RA-STM 
assumes that cracks will develop in the direction parallel 
to the principal compressive stresses in concrete elements, 
and the cracks will "rotate" to follow the principal 
stresses over the entire loading history. In contrast, the 
FA-STM assumes that cracks will develop along the 
direction of principal compressive stresses at initial 
cracking, and the cracks will be "fixed" at this angle 
thereafter. 

The advantage of FA-STM over RA-STM was that 
FA-STM was capable to take into account the concrete 
contribution, induced by the shear stresses along the 
cracks. Nevertheless, the FA-STM models average 
responses, without considering the specific contributions 
of the individual mechanical effects. In this paper, the 
FA-STM is modified by introducing the aggregate 
interlock and dowel action, which are used to control the 
average shear capability of concrete with cracks. In this 
paper, the incorporation of the set of modified FA-STM 
constitutive laws in to the computer code FEAPpv is 
described, and the comparison of the prediction of failure 
load of reinforced high strength concrete with the 
experimental results is given. 



Equilibrium equations 

Assuming that the steel bars can resist only axial 
stresses, then the superposition of concrete stresses and 
steel stresses as shown in Fig. 1 . 
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In the fixed-angle model, the x'-y' coordinate system 
as shown in Fig. 2. is defined. In this coordinate system, 
x and y' are the principal axes of stresses in concrete at 
initial cracking. Angle § is the fixed angle between x and 

x' axes. 
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Fig. 2 Definition of coordinate systems and fixed angle. 

The transformation of stresses in concrete from the x'-y' 
to the x-y coordinate system is given as follows: 
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where, c = cos(-(|)) and s = sin (-(])) . The stresses 
G ct ,,G cv , and T , , are the stresses in concrete in x'-y' 
coordinate system. 

After introducing Eq. 2 into Eq. 1, the final expression 
for equilibrium condition for reinforced concrete can be 
obtained as: 
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Compatibility Equations 

Assuming that no slipping occurs between concrete 
and steel bars, the transformation of the average strains in 
reinforced concrete from the x-y to the x'-y' coordinate 
system is given as follows: 
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where, c = cos ((]) ) and s = sin ((]) ) . The strains e t , , 8 , 
and y,. are the strains in the x'-y' coordinate system. 

Cracking Criterion 

The constitutive relationships of concrete must be 
guided by an interactive cracking criterion for concrete. 
A cracking criterion as shown in Fig. 3 is given as 
follows: 
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where G lc and G 2c is the principal stress in concrete. 
The uniaxial tensile strength f is defined as 
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Fig. 3 Cracking surface. 



As for the cracking envelope under biaxial stress, 
derived for the tension-compression domains and the 
Aoyagi-Yamada model 7 for the domains of tension- 
tension together with the Kupfer's model 7 for the 
compression-compression domains are adopted. 

Constitutive Relationships Of Concrete Before 
Initial Cracking 

Before initial cracking, assuming that the principal 
direction of stress in concrete is coincide with the 
principal direction of strain, the constitutive relationships 
of concrete are given as follows: 



G lc = E c £ 1 , tension 



°" 2c = // 
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where E c is the modulus of elasticity of concrete, 8i and 82 
are the strains of concrete in the ascending and 
descending branches, respectively, and e t ' is the strain of 
concrete at peak stress. 

Constitutive Relationships of Concrete After 
Initial Cracking 

After initial cracking, the constitutive relationships of 
concrete are established in the x'-y' coordinate. 



© ITS JOURNAL OF CIVIL ENGINEERING / Vol. 29 No. 1/ May 2009 



11 



Concrete in compression-tension 

After initial cracking, the stress and strain softening 
occurs in concrete in compression-tension domains. 
The average stress-strain curve of concrete in 
compression 5,6 ' 9 as shown in Fig. 4 is expressed as: 



<V=C/ C ' 



4/C -1 
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(7b) 



where f c is the cylinder compression strength of 
concrete; e c is the concrete strain at maximum 
compressive stress; and £ is the softening coefficient. In 
the descending portion of the concrete stress-strain curve 
the lowest stresses value was taken as 0.2^ f c to avoid 
the potential numerical problem in calculation. 
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Fig. 4. Softened compressive stress-strain curve of 
concrete. 

In Eq. (7a) and Eq. (7b), the stress-softened 
coefficient and the strain-softened coefficient are the 
same value of £ , which can be expressed conservatively 

as: 5 ' 6 
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where £ x , is the tensile strain at x'-direction; p st , p are 

the reinforcement ratios in the x and y directions, 
respectively; f sxY , f syY are the yield stress of steel in the x 

and y directions, respectively; and O x ,a y are the applied 
stresses in the x and y directions, respectively. The 
parameter r|' is less than unity. 



Concrete in tension 

The average stress-strain curve of concrete in tension 1 
as shown in Fig. 5, can be expressed as: 
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where E c is the modulus of elasticity of concrete; /, is 
the uniaxial tensile strength as defined in Eq. (5); e t0 is 

f 

the cracking strain of concrete, which equals to — 

E c 

When concrete is cracked by tension, the tensile stress 
can drop fast. Especially in plain concrete, the tensile 
stress can change much faster than reinforced concrete 
because the crack in concrete can grow without the 
restriction by the reinforcing bars. 

Concrete in shear 

A rational shear modulus has been derived by Zhu. 1 1 
The average stress-strain relationship of concrete in shear 
can be expressed as: 
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Transfer of shear forces across the crack surface in 
reinforced concrete member may result in a large sliding 
deformation and final failure by shear. Nevertheless, the 
model is "smeared". It models average responses, without 
considering the specific contributions of the individual 
mechanical effects, such as the aggregate interlock and 
dowel action at the crack location. 
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Average tensile stress-strain curve of concrete. 



Aggregate Interlock 

Sliding shear deformations at the cracks in reinforced 
concrete structures are resisted by the aggregate interlock 
resulting from roughness of the crack faces. Aggregate 
interlock force is a major factor attributed by concrete in 
the analysis of shear behaviors in reinforced concrete 
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beams. This is mainly due to their significant contribution 
in shear transfer along crack interface of reinforced 
concrete. It is known that at the cracks of a reinforced 
concrete element subjected to in-plane loads, only steel 
forces are initially developed. Aggregate interlock forces 
are produced later to help steel forces sustain the external 
loads. 

Many experimental investigations have been carried 
out for aggregate interlock. Based on both test results and 
the physical concepts of aggregate interlock behavior, 
Bazant and Gambarova 12 developed a more 
comprehensive model. Unlike previous studies, their 
model accounts for the facts discovered in the test: (1) at 
zero sliding-shear deformation, no normal stress exists 
across an open crack; (2) at a constant value of sliding- 
shear deformation, the shear stress tend to decrease as the 
crack opening increases; (3) the shear stresses tend to 
increase with increasing sliding-shear defor-mation at a 
constant crack width. Their ex-pressions for aggregate 
interlock are shown here: 



x a =0.245/; 



0.01D 
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10//; + 2.44(l-16.31//;)|r|' 
1 + 2.44(1- 16.31 1 f' c )r 4 
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where, f' c is the uniaxial compressive strength of 
concrete; 8 n is the crack opening; 8, is the crack slip; 
r = 8, / S n ; and D is the maximum aggregate size (mm). 

In concrete containing many parallel cracks, the 
deformations due to the cracks may be considered 
continuously distributed or smeared. The opening and slip 
of the crack can be expressed in the form of average 
strains due to the smeared cracks (see Fig. 6): 
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where, e x , is the average strain at the crack normal 
direction; and y x . . is the average shear strain at the crack 
direction. 




Fig. 6. Meaning of Symbols in Eq. (11). 

Here l s is the average spacing of the diagonal cracks 
which is calculated from the CEB-FIP Code: 
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where, (]) is the initial cracking angle between the x- 
direction reinforcement and the direction of the initial 
cracking principal tensile stress (see Fig. 2); and 

are the crack spacing in the two orthogonal directions. 
The values of l sx , I are estimated as: 
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The meanings of the symbols in above equations are 
shown in Fig. 7. The symbol s is the spacing of web 

reinforcement, and p sx , p are the reinforcement ratios 

of x- and y-directions. k is a characteristic factor of 
reinforcement and taken as 0.4 for deformed bars or 0.8 
for plain bars or bonded strands. 
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Fig. 7. Meaning of Symbols in Eq. 14a and 14b. 

Based on Bazant's rough cracks model, by optimizing 
the fits of Paulay and Loeber's test data 13 and of Daschner 
and Kupfer's test data, 14 and by assuming for crack 
roughness the aggregate grading suggested by Fuller's 
curve, the following equation was formulated: 15 



(15a) 




x a =0.25/; 

9.8 ///+ 2.44(1- 16 /f')\rf 
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where, r, 8, and 8 n have the same meanings as Bazant's 
equation. However, by comparison with the experimental 
data, it is found that the Eq. 15a may overestimate the 
aggregate interlock of high strength concrete. In the 
proposed model, the equation of aggregate interlock is 
modified to; 
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Dowel Action at the Crack Location 

Dowel action is one of the main mechanisms of load 
transfer along reinforced concrete interfaces. There are 
two possible failure modes of the dowel mechanism. One 
is due to yielding of the dowel bar and the crushing under 
the dowel, the other one is due to the concrete splitting. 
Generally, concrete cover c is the main parameter upon 
which the mode of the failure of the dowel mechanism 
depends. 

The previous experimental investigations 16 have 
shown that when c is greater than 6 to 7 time the bar 
diameter, failure is controlled by the crushing of concrete 
and the yielding of the bar. For smaller concrete cover, 
the failure mechanism is governed by splitting of the 
concrete, and the splitting cracks can be opened at the 
bottom or at the side faces of the section. In a reinforced 
concrete element or a beam, the two directions 
reinforcement is usually arranged along the brim of the 
cross section, thus these bars have smaller side covers, 
and dowel failure occurring in these structural members 
belongs to the concrete splitting. 

In previous work, 17 ' 2 the existence of the shear forces 
due to dowel action was neglected, and this 
underestimated the role of the steel. In this present 
research, the shear contributions from dowel action are 
considered until the final collapse of the element. 

Generally, the dowel action and aggregate interlock 
have some common characteristics. They both occur at 
the crack location in order to replace the lost stress in the 
steel and concrete. Both are transmitted to the bulk of the 
concrete, and therefore, they are in fact concrete forces. 

For the prediction of the dowel action when failure is 
due to the concrete splitting, only the empirical equations 
are available. Eq. 16 was developed by Baumann and 
Ruch 17 . 



D.=lMb n d b l[f: inN,mm 



(16) 



where, b n is the net width of beam and d b is the diameter 
of the dowel bar (see Fig. 8). 

The average stresses should be limited by the stresses 
at the crack. Therefore, the average shear stress should be 
controlled by the local shear stress contributed by the 
aggregate interlock and dowel action. 
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Fig. 8. Meaning of Symbols in Eq. 16. 
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Fig. 9 Average stress-strain curve of steel bars embedded 
in concrete. 

Constitutive relationship of steel bar 

The stress-strain curve of steel can be modeled by two 
straight lines, 2 ' 4 ' 10 as shown in Fig. 9. The bilinear model 
is given as follows: 



° s = fy 



a = E e , e <e 

s s s ' s n 

(0.91-2fi) + (0.02 + 0.25fl)- 



(18a) 
(18b) 



where B is a parameter defined as (l / p ) (f cr I f Y ) , and 
e„ is the average yield strain of mild steel bars embedded 
in concrete at the beginning of yielding, taken as 
e r (0.93-2B). 



ft xY | < ft u | + (| D x sin (p | + 1 D y cos (p |) / bd v (17) 

where, D x , D y are the dowel action of all of y- of x- 

direction bars, respectively, b is the width of a 

rectangular beam, d v is the depth of effective shear 
element, and § is the initial cracking angle. 



Test Program 

Eighteen simply supported reinforced high-strength 
concrete beams were to be tested under two-point loading. 
It is already known that the thickness of the beams has no 
contribution to size effect, so all test beams were chosen 
to have the same width of 185 mm. The test beams were 
divided into two series with different a/d ratio of 2 and 
3.5. 
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Table 1. Details of 18 experimental beams 



Groups Beam Main tension Vertical web 





Notations 




Beam dimension 




Concrete strength 


Steel strength 


reinforcement 


reinforcement 






b 


h 


d 


a 


1 


fc' Batch 


fys 


fyv 


As 


P 


Av 


s 


pv 






mm 


(mm) 


(mm) 


(mm) 


(mm) 


MPa 




MPa 


MPa 


mm2 


% 


mm2 


mm 


% 


Series 1 




B-2-200 


185 


250 


200 


400 


1300 


89 


1 


512 


0 


1480 


3.98 


- 


- 


0 


1 


B-2-400 


185 


475 


400 


800 


2600 


95 


1 


512 


0 


2960 


3.98 


- 


- 


0 




B -2-700 


185 


825 


700 


1400 


4550 


92 


3 


512 


0 


5180 


3.98 


- 


- 


0 




V-2-200 


185 


250 


200 


400 


1300 


75 


4 


512 


518 


1480 


3.98 


57 


91 


0.34 


2 


V-2-400 


185 


475 


400 


800 


2600 


95 


1 


512 


518 


2960 


3.98 


157 


250 


0.34 




V-2-700 


185 


825 


700 


1400 


4550 


98 


3 


512 


518 


5180 


3.98 


157 


250 


0.34 




VV- 2-200 

V V Z. iUU 


185 


250 


200 


400 


1300 


75 


4 


512 


518 


1480 


3.98 


157 


150 


0.57 


3 


VV- 2-400 


185 


475 


400 


800 


2600 


95 


1 


512 


518 


2960 


3.98 


157 


150 


0.57 




V V-2-700 


185 


825 


700 


1400 


4550 


98 
Series 2 


3 


512 


518 


5180 


3.98 


157 


150 


0.57 




B-3.5-200 


185 


250 


200 


700 


1900 


85 


2 


512 


0 


1480 


3.98 






0 


1 


B-3.5-400 


185 


475 


400 


1400 


3800 


85 


2 


512 


0 


2960 


3.98 






0 




B-3.5-700 


185 


825 


700 


2450 


6650 


102 


1 


512 


0 


5180 


3.98 






0 




V-3.5-200 


185 


250 


200 


700 


1900 


85 


2 


512 


518 


1480 


3.98 


57 


91 


0.34 


2 


V-3.5-400 


185 


475 


400 


1400 


3800 


85 


2 


512 


518 


2960 


3.98 


157 


250 


0.34 




V-3.5-700 


185 


825 


700 


2450 


6650 


102 


1 


512 


518 


5180 


3.98 


157 


250 


0.34 




VV-3.5-200 


185 


250 


200 


700 


1900 


85 


2 


512 


518 


1480 


3.98 


157 


150 


0.57 


3 


VV-3.5-400 


185 


475 


400 


1400 


3800 


85 


2 


512 


518 


2960 


3.98 


157 


150 


0.57 




VV-3.5-700 


185 


825 


700 


2450 


6650 


104 


1 


512 


518 


5180 


3.98 


157 


150 


0.57 



Notations: B-2-200 = "B" means beams without web reinforcement; 2 means a/d ratio of 2; 200 is effective depth,; "V" 
means beams with web reinforcement percentage of 0.34, "VV" means beams with web reinforcement 
percentage of 0.57, etc. b — width of the beam, h — overall depth of beam, d — effective depth of beam, a — 
shear span, I — overall length of beam, p — main tension steel ratio, p v — web reinforcement ratio. 



Table 2. Cracking loads, ultimate Loads and failure modes of 18 concrete beams 



No. 


Beam 


Cracking Stress (MPa) 


Ultimate 


Modes 


Ratios of Loads 


of 




Flexural 


Diagonal 


Stress (Mpa) 


of 






Beams 


Notations 


vf 


vd 


vu 


Failure 


vf /vu 


vd/vu 


1 


B-2-200 


0.68 


203 


6.28 


DS 


0.10 


0.32 


2 


V-2-200 


0.41 


2.03 


7.57 


DS 


0.05 


0.27 


3 


VV-2-200 


0.54 


2.7 


8.23 


CC 


0.07 


0.33 


4 


B-3.5-200 


0.41 


1.89 


1.89 


DT 


0.22 


1.00 


5 


V-3.5-200 


0.27 


1.62 


4.32 


DT 


0.06 


0.38 


6 


W-3.5-200 


0.41 


1.89 


4.69 


CC 


0.09 


0.40 


7 


B-2400 


0.74 


1.89 


6.49 


DS 


0.11 


0.29 


8 


V-2-400 


0.95 


2.03 


8.11 


DS 


0.12 


0.25 


9 


VV-2-400 


0.95 


2.36 


8.23 


CC 


0.12 


0.29 


10 


B-3.5-400 


0.27 


1.82 


1.82 


DT 


0.15 


1.00 


11 


V-3.5400 


0.41 


1.22 


5.03 


DT 


0.08 


0.24 


12 


W-3.5400 


0.41 


1.62 


5.07 


CC 


0.08 


0.32 


13 


B-2-700 


0.77 


1.54 


4.44 


DS 


0.17 


0.35 


14 


V-2-700 


0.62 


1.54 


6.22 


DS 


0.10 


0.25 


15 


VV-2-700 


0.62 


1.54 


7.61 


sc 


0.08 


0.20 


16 


B-3.5-700 


0.39 


1.54 


1.54 


DT 


0.25 


1.00 


17 


V-3.5-700 


0.39 


1.54 


4.4 


DT 


0.09 


0.35 


18 


W-3.5-700 


0.46 


1.54 


4.79 


CC 


0.10 


0.32 



Notations: DS - diagonal splitting failure; DT - diagonal-tension failure; FC - flexural failure in 
compression zone; SC - shear compression 



In each series, three groups of beams were cast 
according to the different percentages of web 
preinforcement. In each group the effective depths of the 
specimens were varied from 200 to 700 mm, whereas all 
other variables remained constant. The specimen 
dimensions and other parameters are given in Table 1 . 



The beams were geometrically scaled appropriately in 
all respects (see Fig. 10) except the maximum aggregate 
size which is maintained at 20 mm. Previous investigation 
by Walraven and Lehwalter 19 indicated that there was no 
significant influence of the maximum aggregate size on 
the size effect in shear behavior. All of the specimens 
were designed to fail in shear as far as possible. 
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Fig. 11. Typical test set-up. 
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Notes: Three beams for each figure, 
one beam without stirrups, 
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Fig. 10. Details of beam specimens. 

Materials 

The compressive strength of the concrete was 
designed to be about 100 MPa for all the beam specimens. 
The slump of the mix was about 200 mm. Two 150 x 150 
x 150 mm concrete cubes and two 150 x 300 mm 
concrete cylinders were cast and tested together with each 
specimen to obtain the compressive strength of concrete, 
f cu and f' c . Since the beams were cast in different 
batches, regression analyses were carried out for concrete 
compressive strength in each batch. Different average 
concrete strength of each beam specimen was then 
obtained from the results of these regression curves 
according to the age of specimen at the time of testing. 
The results are given in Table 1 . 

There are two types of reinforcement bars used in the 
beam specimens: deformed high tensile steel (T bar) and 
plain round mild steel (R bar). The average yield stresses 
of the steels used in the specimens are tabulated in Table 
1. Note that T10 indicates a T bar of 10 mm diameter, etc. 

Test set-up 

The typical set-up for specimens is shown in Fig. 11. 
The test specimens were loaded by two 1000 kN actuators. 
Before testing, each specimen was white-washed on one 
surface and 150 mm x 150 mm grids were drawn in 
pencil to facilitate crack detection. Deflections of the 
beam specimens were measured using Linear Voltage 
Displacement Transducers (LVDTs) located at the bottom 
central line of the beams at seven locations. 
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COMPARISON BETWEEN THEORETICAL 
ANALYSIS AND EXPERIMENTAL DATA 

To verify the reliability of the proposed method, the 
18 beams are analyzed as shown in Table 2. The 
comparisons between the proposed model and the 
experimental data are shown in Figs. 12-17 that show the 
relationship between the load and midspan deflection. 
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The calculation was terminated when the average 
shear strain in the element reaches 1.0%. It was found that 
the shear strain increases suddenly afterwards, leading to 
failure. In these figures, it is found that the load- 
displacement curves calculated by the model without 
considering the aggregate interlock and dowel action 
almost coincide with the curves calculated by the model 
with considering the aggregate interlock using Eq. 15a. 
Eq. 15a may overestimate the aggregate interlock of high 
strength concrete, especially for the high strength concrete 
deep beams without web reinforcement as shown in Fig. 
15(a), Fig. 16(a) and Fig. 17(a). So Eq. 15a is changed 
into Eq. 15b to reduce the overestimation of aggregate 
interlock. 
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Fig. 17. Load-displacement relations of deep beams with 
(jf=700mm (continued). 



CONCLUSIONS 

In this paper a modified FA-STM for the analysis of 
the behavior of high strength reinforced concrete beams, 
subjected to plane stresses and monotonically loaded up 
to failure, is proposed. The crack angle is kept fixed at the 
initial direction of concrete. The functions for aggregate 
interlock and dowel action are used to control the shear 
capacity of cracked concrete. The reliability of the model 
has been assessed through the comparison with the 
experimental results. The comparison shows the 
aggregate interlock and dowel action can reduce the 
overestimation of the shear capacity of high strength 
reinforced beam, especially the high strength reinforced 
deep beam without web reinforcement. Moreover, the 
model is suitable for being implemented numerical 
procedures due its simplicity. 
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